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1.0 


INTRODUCTION 


During  the  past  three  years,  under  joint  US  Army  Research 
Office  and  The  University  of  Tennessee  funding, personnel  of  the 
Gas  Diagnostics  Research  Division  at  The  University  of  Tennessee 
Space  Institute  have  developed  a  unique  three-component  laser  velo- 
cimeter  for  the  in  situ  measurement  of  particle  and/or  gas  velocities 
in  flow  fields  produced  behind  bodies  at  high  angles  of  attack  and 
in  jet  exhaust  plumes. 

This  report  describes  the  development  of  the  laser  veloci- 
meter  and  its  subsequent  application  to  the  measurement  of  the  velocity 
distribution  and  vortex  structure  in  free  jets  and  in  flows  over 
missiles  at  high  angles  of  attack. 


2.0  TECHNICAL  OBJECTIVES  AND  APPROACH 


2.1  Technical  Objective 

The  overall  technical  objective  of  the  research  program  was 
to  develop  and  utilize  a  three-component  Laser  Doppler  Velocimeter (LPV)  to 
measure  the  velocity  distribution  and  vortex  structure  in  flows  over 
missiles  at  high  angles  of  attack.  Three  dimensional  velocity  measure¬ 
ments  were  needed  to  obtain  a  better  characterization  and  understanding 
of  the  complex  gas  dynamic  phenomena  and  flow  structures  downstream 
of  bodies  at  high  angles  of  attack. 

In  addition,  the  measurements  would  provide  detailed  quantitative 
data  to  verify  existing  and  future  analytical  models  for  the  design  of 
missile  systems. 


2.2  Technical  Approach 

In  the  first  phase  of  the  research  program  outlined  in  Table  2.1, 
a  2D  laser  velocimeter  system  was  to  be  designed  and  constructed  which 
could  be  adapted  to  provide  a  three-dimensional  capability  during  the 
program.  As  described  later  in  Section  3.1,  a  two-component,  dual- 
scatter,  crossed-beam  laser  velocimter  optical  system  was  developed 
and  used  to  obtain  the  velocity  distributions  in  the  wake  of  a  circular 
cylinder  and  in  the  exhaust  plume  produced  by  a  subsonic  nozzle. 

To  provide  flexibility  in  the  experimental  setup  and  application, 
the  laser  velocimeter  signal  processing  system  was  designed  using  a 
fully-programmable  microprocessor.  This  system  is  described  in 
Section  3.1.2. 


At  the  beginning  of  the  program,  it  was  anticipated  that  the 
2D  laser  velocimeter  system  would  be  calibrated  by  comparing  LDV  data 
with  hot  wire  anemometer  data  obtained  in  the  exhaust  flow  field  pro¬ 
duced  by  a  small  subsonic  nozzle  discharging  into  the  atmosphere. 

Upon  completion  of  the  calibration  phase,  the  2D  LDV  system  would  be 
used  to  measure  the  velocity  distribution  downstream  of  a  circular 
cylinder  placed  perpendicular  to  the  major  flow  direction.  As  described 
in  Section  5.1,  this  was  accomplished  as  planned. 

As  indicated  in  Table  2.1,  the  next  task,  Task  III,  consisted 
of  extending  the  2D  LDV  system  to  a  3D  system.  It  was  envisioned  that 
the  third  velocity  component  would  be  measured  using  a  reference  beam 
technique.  Considerable  time  and  expense  was  devoted  to  the  development 
of  a  HeNe  laser  reference  beam  optical  system  (described  in  Section  3.2) 
which  would  scan  in  conjunction  with  the  scannable  2D  dual  scatter 
optical  system.  A  unique  afocal  scan  system  was  developed  which  enabled 
the  laser  probe  volumes  to  be  coincident  and  to  be  scanned  through  a 
flow  field  of  Interest.  Unfortunately,  the  inadequate  spatial  resolution 
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of  the  third-component,  reference  beam  system  resulted  in  the  modi¬ 
fication  of  the  experimental  setup.  The  final  system  employed  to 
obtain  the  three-component,  three-dimensional  velocity  distribution 
data  utilized  two  two-component  dual-scatter  systems.  The  first 
system  measured  two  velocity  components  (usually  u  and  v)  while  the 
second  system  also  measured  two  velocity  components  (typically  v  and  w) . 
Due  to  cost  limitations,  the  two  systems  utilized  the  same  digital 
microprocessor  systems  which  means  the  four  velocity  component  measure¬ 
ments  were  not  obtained  simultaneously.  While  this  does  not  affect 
the  mean  velocities,  it  does  mean  that  the  turbulent  cross  correlations 
were  not  obtained  at  the  same  time. 

The  experimental  approach  consisted  of  making  four-component 
velocity  measurements  at  a  point  in  the  flow  field  established  about 
a  body  of  revolution  or  in  a  jet  exhaust.  The  laser  probe  volume  was 
then  translated  in  three  orthogonal  directions  to  obtain  the  spatial 
velocity  distribution. 

As  shown  in  Table  2.1,  Tasks  IV  through  VI,  measurements  were 
to  be  performed  in  a  series  of  tests  starting  at  low  subsonic  Mach 
numbers  in  a  VSTOL  subsonic  wind  tunnel,  at  transonic  Mach  numbers 
in  the  UTSI  Transonic  Wind  Tunnel,  and  at  low  supersonic  Mach  numbers 
in  the  NASA  Marshall  7-inc.h  Wind  Tunnel.  Due  to  delays  in  the  develop¬ 
ment  of  the  instrumentation,  three  objectives  were  not  completed. 


Task  I 
Task  II 


Task  III 
Task  IV 


Task  V 


Task  VI 


TABLE  2.1 

PROPOSED  RESEARCH  PROGRAM  TASKS 

Design  aid  Assembly  of  2D  Laser  Velocimeter  System. 

Initial  Testing  and  Calibration  2D  Subsonic  Flows 
Over  Circular  Cylinders. 

Extension  of  LDV  System  to  3D  Capability. 

Measurements  of  Low  Subsonic  Flow  Over  Cone-Cvlinder 
Combination  at  High  Angles  of  Attack. 

Extension  of  Laser  Velocimeter  Measurements  to  Compressible 
Subsonic  and  Transonic  Mach  Number  Regimes. 

Laser  Velocimeter  Measurements  at  Low  Supersonic  Mach  Numbers 


Task  VII 


Comparison  of  Experimental  Data  with  Analytical  Models. 


3.0  DESCRIPTION  OF  LASER  VELOCIMETER  SYSTEM 


The  Gas  Diagnostics  Research  Division  has  in  operation  a 
three-component  laser  velocimeter  System  which  is  based  on  the  use 
of  two,  two-component  Bragg-cell  fringe  type  laser  velocimeter  optical 
systems . 

The  LV  signal  processing  system  was  designed  around  a  unique 
Z80  microprocessor  system  which  is  capable  of  real-time,  on-line  data 
reduction.  The  microprocessor  controls  the  laser  velocimeter  signal 
processors  and  traversal  systems  and  is  equipped  with  software  programs 
which  allow  the  recording  of  various  system  parameters  in  addition  to 
the  measured  data. 

Large  amounts  of  data  in  a  histogram  format  can  be  acquired  in 
a  very  short  time  interval.  The  data  include  velocity  histograms 
covering  four  decades  of  velocity  variation  and  particle  size  histograms. 
The  microprocessor  processes  the  data  to  yield  average  velocity  and 
particle  size,  turbulent  intensity  components  and  the  kurtosis  of 
the  signal. 

The  two-component  laser  velocimeter  optical  and  electronic 
systems  are  described  in  Section  3.1.  The  use  of  two,  two-component 
laser  velocimeter  systems  to  obtain  three-dimensional  velocity  distri¬ 
butions  are  described  in  Section  3.2. 


3.1  Two-Component  Laser  Ve loc imeter  System 

Measurement  of  two  true  components  of  velocity  was  provided  bv 
a  double  Bragg  cell  laser  velocimeter  (Ref.  1).  The  optical  configuration 
is  shown  schematically  in  Fig.  3.1.  A  photograph  of  the  svstem  is 
given  in  Fig.  3.2.  A  block  diagram  of  the  electronic  instrumentation 
is  shown  in  Fig.  3.3.  A  noteworthy  feature  of  this  tvpe  of  velocimeter 
is  that  both  components  of  velocity  are  obtained  from  a  single  photo¬ 
multiplier  tube.  Optical  simplicity  is  obtained  at  the  cost  of  in¬ 
creased  electronic  complexity,  but  this  tradeoff  is  desireable  in 
most  practical  applications  of  laser  velocimetrv  because  electronic 
components  withstand  harsh  environments  better  than  optical  components. 


3.1.1  Two-Component  La s e r  Veloc imet er  Optical  Sv s t em 

As  shown  in  Fig.  3.1,  the  TEMor)  beam  from  a  2-watt  Argon  Ton 
laser  is  passed  through  a  double  Bragg  cell  (i.e.  acousto-optic 
modulator)  whose  function  is  to  simultaneously  diffract  and  frequency 
shift  75  percent  of  the  laser  beam  into  three  beams  which  are  angularly 
separated  with  respect  to  each  other  and  the  remaining  25  percent  of 
the  original  beam.  The  Bragg  cell  is  a  small  water  tank  in  which  two 
quartz  crystal  transducers  are  mounted  at  right  angles  with  respect 
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Block  diagram  of  laser  velocimeter  signal  processing  system 


to  each  other.  A  phase-lock-loop  (PLL)  oscillator  drives  a  power 
amplifier  which  in  turn  drives  a  quartz  transducer.  The  frequency 
of  each  PLL  is  tunable  in  order  that  the  transducer  be  operated  at 
its  maximum  efficiency.  The  output  power  of  each  power  amplifier 
is  variable,  and  it  is  this  control  which  allows  one  to  obtain  four 
beams  of  equal  intensity  from  the  double  Bragg  cell. 

As  shown  in  Fig.  3.1,  lens  LI  focuses  the  four  beams  to  a  common 
region  in  space  (the  LV  probe  volume)  producing  two  orthogonal  sets 
of  moving,  planar  interference  fringes.  Light  scattered  from  these 
fringes  by  a  moving  particle  is  collected  by  lens  L2  and  directed  to 
the  photomultiplier  tube  (PMT) . 


3.1.2  Microprocessor  Laser  Velocimeter  Signal  Processing  Electronics 

As  shown  in  Fig.  3.2,  there  are  three  main  components  of  the 
laser  velocimeter  signal  processing  electronics:  the  double  Bragg  cell 
electronics,  the  burst  counter  laser  velocimeter  processor,  and  the 
stepping  motor  controller  electronics.  A  photograph  of  the  Bragg  Cell 
and  Scanner  Electronics  svpttems  is  given  in  Fig.  3.4  while  a  photograph 
of  the  Microprocessor  Lata  Acquisition  System  is  given  in  Fig.  3.5. 

One  of  the  Bragg  cell  transducers  is  operated  at  a  frequency 
of  15  MHz  while  the  other  is  operated  at  45  MHz.  The  power  spectrum 
of  the  PMT  signal  produced  by  a  particle  at  rest  would  show'  two  strong 
spectral  components,  one  at  15  MHz  and  the  other  at  45  MHz  (Ref.  2). 
Motion  of  the  particle  will  shift  these  two  spectral  components  of  the 
signal  away  from  their  center  frequencies.  As  long  as  the  particle's 
velocity  is  not  so  large  that  the  two  spectral  components  of  the  signal 
overlap,  the  two  signal  frequencies  can  be  electronical lv  separated 
and  recorded  individually.  Figure  3.2  presents  a  block  diagram  showi'n.^ 
how  the  two  signal  frequencies  are  separated  and  recorded. 

The  signal  splitter  divides  the  total  signal  amplitude  into  two 
equal  parts,  each  of  which  is  passed  to  a  signal  processing  channel 
consisting  of  a  double  balanced  mixer,  PLL  local  oscillator,  low  pass 
filter,  and  burst  counter  LV  signal  processor.  The  signal  channel 
responsible  for  processing  signals  produced  by  the  15.0  MHz  Bragg 
transducer  is  called  the  15.0  MHz  channel,  and  the  other  signal  channel 
is  called  the  45.0  MHz  channel.  The  LO  (local  oscillator)  port  of  each 
mixer  is  driven  by  a  PLL  whose  frequency  is  adjustable.  For  each 
signal  frequency  which  appears  in  the  mixer's  input  (its  radio  frequence, 
or  RF,  port)  the  mixer's  output  (its  intermediate  frequence,  or  IF.  port) 
contains  the  sum  of  this  frequency  and  the  local  oscillator  frequence, 
and  the  difference  between  this  signal  frequence  and  the  local  oscillator 
frequency.  However,  the  mixer  output  is  passed  through  a  low  pass  filter 
which  lets  through  only  those  difference  frequencies  which  fall  within 
the  pass-band  of  the  filter.  Ray,  for  example,  that  motion  of  a  particle 
through  the  interference  fringes  produces  signal  frequencies  of  16.0  and 
49.0  MHz,  and  that  the  local  oscillator  for  the  15.0  MHz  channel  is  set 
to  17.0  MHz  while  the  45.0  channel  MHz  LO  frequency  is  50.0  MHz.  Assuming 
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Chat  the  low  pass  filter  for  the  15.0  MHz  channel  is  set  to  pass  all 
frequencies  below  2.0  MHz  while  the  low  pass  filter  for  the  45.0  MHz 
channel  Is  set  to  pass  frequencies  below  5.0  MHz,  then  the  output 
frequency  of  the  15.0  MHz  channel  will  he  1.0  MHz  whereas,  the  output 
of  the  45.0  MHz  channel  will  be  4.0  MHz. 

The  signal  frequency  produced  by  a  particle  at  rest,  f0,  is 
the  difference  between  the  frequency  of  the  local  oscillator,  fpo, 
and  the  Bragg  excitation  frequency,  fg, 

fo  “  fL0  “  ‘B 

A  component  of  velocity  for  a  moving  particle  Is  found  from 
the  difference  between  fQ  and  the  frequency  of  the  output  of  the  low 
pass  filter,  f.  The  latter  is  measured  by  a  burst  counter  LV  signal 
processor  (actually,  the  burst  processor  measures  T,  the  signal  period 
averaged  over  some  integer  number  of  periods).  Velocity  is  computed 
from  the  equation 


v  =  6 (fa  -  f ) 

in  which  6  is  the  fringe  period  (i.e.  distance  from  one  interference 
fringe  to  the  next).  The  magnitude  of  a  component  of  velocity  is 
determined  from  the  absolute  value  of  fp  -  f,  and  the  sign  of  the 
component  of  velocity  is  determined  from  the  sign  of  f 0  -  f . 

A  "simultaneity  checker"  rejects  signal  bursts  that  do  not 
at  any  point  in  time  simultaneously  appear  in  both  signal  channels. 
Thus,  the  magnitudes  and  signs  of  two  of  the  three  components  of 
velocity  are  measured. 

The  burst  counter  TV  signal  processors  are  constructed  in  the 
form  of  interface  cards  which  mount  in  a  7.80  based  computer.  Each 
processor  contains  two  fringe  counters  which  are  started  simultaneously 
but  are  otherwise  independent.  One  fringe  counter  converts  an  integer 
number  of  signal  periods  (called  the  short  fringe  count)  into  a  single 
pulse  lasting  that  many  signal  periods.  The  second  fringe  counter 
converts  a  larger  number  of  signal  periods  (called  the  long  fringe 
count)  into  a  pulse  lasting  this  many  signal  periods.  The  long  and 
short  fringe  count  pulses  gate  the  output  of  a  100.0  MHz  oscillator 
to  the  inputs  of  the  long  and  short  clock  counters,  respectively. 

If  N  is  a  fringe  count  and  n  is  the  number  of  clock  periods  counted 
while  the  fringe  count  was  being  acquired,  then  the  signal  period  is 
given  by 


wherein  our  case  the  clock  frequency,  f  ,  ,  is  100.0  MHz.  Such  a  computation 
can  be  performed  for  both  the  long  and  snort  fringe  counting  processes, 
and  they  should,  of  course,  give  substantially  the  same  result  for  the 
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signal  period.  In  practice,  the  two  computations  do  not  always  agree 
and  the  maximum  allowed  disparity  between  the  two  is  called  the 
"aper iodicity  limit."  Data  for  which  the  aperiodicity  between  the 
long  and  short  fringe  counts  exceeds  this  limit  are  rejected.  A 
typical  value  of  the  aperiodicity  limit  is  3.0  percent. 

An  LV  signal  must  initially  rise  in  amplitude  above  the  system 
electronic  noise  and  it  will  eventually  disappear  back  into  the  noise. 
Hence,  the  initial  and  final  signal  periods  of  the  signal  burst  are 
unavoidably  of  poor  signal-to-noise  ratio.  A  third  fringe  counter, 
called  the  precounter,  preceeds  the  long  and  short  fringe  counters. 

The  function  of  the  precounter  is  to  prevent  the  long  and  short  fringe 
counters  from  counting  the  first  few  signal  periods.  The  precount  is 
settable  from  1  to  15. 

As  a  result  of  various  signal  conditions,  the  signal  mav  drop 
to  zero  before  the  fringe  counters  have  completed  their  prescribed 
counts.  Unless  a  means  of  detecting  this  condition  is  designed  into 
the  LV  signal  processor,  the  instrument  will  spend  almost  all  of  its 
time  in  this  "signal  dropout"  condition.  The  signal  processor  contains 
a  zero  crossing  detector  (ZCD)  whose  function  is  to  convert  the  analog 
input  to  a  burst  of  digital  pulses  having  the  same  time  period  as  the 
signal  period.  While  processing  a  signal  burst  the  ZCD  output  shou1d, 
on  the  average,  be  logically  TRUE  about  the  same  percentage  of  time 
that  it  is  logically  FALSE.  Our  dropout  detector  consists  of  a  circuit 
which  measures  the  average  ratio  of  the  ZCD  output  TRUE  time  to  its 
FALSE  time  for  each  signal  burst.  If  the  FALSE  time  exceeds  the  TRUE 
time  by  a  factor  of  two,  the  processor  is  immediately  reset  in  pre¬ 
paration  for  the  next  burst.  Recovery  from  a  brief  burst  of  noise 
requires  only  a  few  microseconds. 

An  S-100  bus  Z80  computer  controlls  essentially  all  svstem 
operating  parameters.  The  noteable  exceptions  are  the  PMT  high 
voltage,  the  laser  power,  the  frequencies  of  the  four  PLL  oscillators 
(two  of  which  drive  the  double  Bragg  cell  while  the  other  two  serve 
as  local  oscillators  for  the  mixers),  and  the  v  (vertical)  and  the 
z  (downstream)  spatial  positions  of  the  LV  probe  volume.  Spatial 
scanning  in  the  x  (transvers)  direction  was  provided  under  computer 
control  ot  two  stepping  motor  driven  translation  stages,  one  of  which 
carried  lens  LI  (nee  Fig.  3.1)  while  the  second  carried  lens  1.2  and  the 
PMT  together.  A  complete  transverse  profile  could  be  recorded  without 
operator  intervention. 

A  timer  having  a  resolution  of  from  1.0  second  to  1.0  microseconds 
is  read  at  the  start  of  each  signal  hurst.  As  a  result,  the  data  obtained 
are  records  of  two  components  of  velocity  versus  time.  These  data  are 
stored  on  8  inch  floppy  diskettes. 
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3_._2  Three  - Comp  on  en  t  Lap  e_r_  Velocl  meter  Optical  Systeri 

As  discussed  earlier,  it  had  been  expected  that  the  UTST 
3D  laser  velocimeter  optical  system  (shown  schematically  In  Fig.  3.b), 
which  employs  an  Afocal  Fean  System, would  be  fabricated  and  utilized 
in  this  research  program.  As  can  be  seen,  the  system  utilized  a  Bragg 
cell  optical  arrangement  to  obtain  the  2D  velocitv  vector  and  a  low 
power  HeNe  laser  Doppler  reference  beam  svstem  to  obtain  the  third 
component.  This  svstem  had  the  advantage  of  operating  total lv  in 
the  backscatter  mode  and  could  he  operated  to  scan  the  laser  probe 
volume  across  the  flow  of  interests. 

Unfortunately,  the  spatial  resolution  for  the  third  component 
proved  to  be  inadequate  for  the  flow  phenomena  associated  with  subscale 
models  in  the  wind  tunnel.  The  I'TST  Afocal  Optical  Fean  Svstem  would 
be  ideal  for  full-scale  operations. 

In  order  to  achieve  a  true  three-component  svstem,  due  to  cost 
limitations,  it  became  necessarv  to  develop  two  two-c omponent  optical 
systems  which  would  utilize  the  same  microprocessor  electronics  systems. 
A  schematic  diagram  of  the  system  s  shown  in  Fig.  3.7.  As  shown  in 
the  diagram,  a  moveable  mirror  was  used  to  translate  the  four  laser 
beams  from  the  Bragg  cell  optics  in  two  directions  relative  to  the 
flow.  With  the  translatable  mirror  in  place,  the  laser  velocimeter 
system  provided  a  measurement  of  the  u  and  v  velocitv  components. 

Without  the  translatable  mirror  in  place,  the  v  and  w  velocitv  com¬ 
ponents  were  measured.  The  same  signal  detector  and  electronics  were 
employed  at  all  times,  so  that  the  velocitv  measurements  were  not 
simultaneous . 
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EXPERIMENTAL  FACILITIES  UTILIZED  IN  MEASUREMENT  PROGRAM 


Due  to  a  lack  of  available  wind  tunnels,  the  measurements 
described  in  this  report  were  performed  in  a  Free  .let  Facility  at 
UTSI.  The  subsonic  free  jet  facility  is  described  in  Section  4.1 
while  the  supersonic  free  jet  facility  is  described  in  Section  4.2 


4.1  UTSI  Subsonic  Free  Jet  Wind  Tunnel  Facility 

A  schematic  diagram  of  the  UTSI  Subsonic  Free  Jet  Wind  Tunnel 
is  given  in  Fig.  4.1.  The  Free  Jet  Facility  utilizes  the  UTSI 
Transonic  Wind  Tunnel  high  pressure  air  storage  capacity  (750  cubic 
feet  at  3200  psia).  The  maximum  flow  rate  is  190  lbm/sec  which  is 
sufficient  to  permit  continuous  operation  of  the  subsonic  free  jet 
nozzle  used  in  the  test  described  in  this  report. 

As  shown  in  Fig.  4.1,  air  at  a  few  hundred  pounds  pressure 
enters  a  large  plenum  chamber  where  it  passes  through  a  solid  particle 
filter  and  baffle  system  before  it  exhausts  into  the  ambient  atmosphere 
through  a  subsonic  (normally  M  =  0.2)  nozzle. 

In  order  to  provide  submicron  size  particles  which  closely 
follow  the  gas  streamlines,  the  flow  is  seeded  with  aluminum  oxide 
(AI2O3)  particles  of  known  submicron  diameter.  The  seed  particles 
are  added  only  as  required  to  increase  the  data  rate.  For  many  tests 
the  submicron  particles  remaining  after  the  filtering  process  wore 
present  in  sufficient  numbers  so  that  no  additional  seed  particles 
were  required. 

A  typical  jet  free  stream  velocity  profile  at  the  nozzle  exit 
plane  is  presented  in  Fig.  4.2  for  a  let  exit  velocity  near  65  m/sec 
which  corresponds  to  a  Mach  number  of  approximately  0.2.  Additional 
velocity  profiles  are  described  in  Section  5.1. 

In  this  type  of  wind  tunnel,  the  model  is  placed  in  the  so-called 
potential  core  region  of  the  free  jet  exhaust.  For  the  exit  Mach  number 
of  0.2,  the  length  of  the  potential  core  is  approximately  20  cm  long 
with  widths  of  30mm  at  the  nozzle  exit  plane  and  approximately  10mm  at 
18mm  from  the  nozzle  exit.  The  size  of  this  "test  section"  greatly 
restricted  the  size  of  the  models  tested  and  prevented  data  collection 
in  the  model  boundary  layer. 
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4.2  UTST  Supersonic  Jet  Exhaust  Facility 


Some  measurements  in  the  exhaust  flow  field  from  a  Mach  2 
axisvmmetric  nozzle  were  obtained  using  the  IITS1  Laser  Velocimeter 
System.  The  schematic  diagram  of  the  Supersonic  Jet  Exhaust  Facility 
is  the  same  as  that  presented  in  Fig.  4.1.  A  photograph  of  the  Mach 
2  Free  Jet  Facility  is  given  as  Fig.  4.3. 

The  jet  nozzle  throat  diameter  is  1.0  inch  while  the  exit 
diameter  is  1.3  inch.  The  exit  pressure  was  maintained  at  25  psia 
which  produced  a  mass  flow  rate  of  3.5  Ib/sec.  The  jet  could  be 
operated  at  this  flow  rate  for  up  to  30  minutes. 

Measurements  of  the  flow  field  produced  in  this  facility  are 
described  in  Section  5.4. 
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Fjg ,  4.1  Schematic  diagram  of  I’TST  Subsonic  Free 
Jet  Wind  Tunnel . 


5.0  DISCUSSION  OF  RESULTS 


5.1  Subsonic  Jet  Exhaust  Plume  Measurements 


The  UTSI  Laser  Velocimeter  System  Operating  in  the  2-D  mode  was 
utilized  to  obtain  the  velocity  distribution  in  the  exhaust  from  a  subsonic 
nozzle  having  an  exit  diameter  of  30  mm.  The  nozzle  exhausted  into  an  ambient 
environment.  A  high  pressure  air  supply  was  used  to  provide  a  continuous  flow 
of  air.  The  air  was  filtered  to  remove  particles  having  diameters  greater  than 
one  micron.  The  remaining  submicron  size  particles  in  the  flow  were  present 
in  sufficient  numbers  for  the  measurements  reported  here.  It  was  assumed  that 
the  submicron  size  particles  followed  the  gas  streamlines  and  moved  at  the 
gas  velocity  except  in  regions  of  locally  high  velocity  gradients. 

The  velocity  data  reported  here  were  obtained  by  programming  the 
laser  velocimeter  microprocessor  to  scan  the  laser  probe  volume  across  the 
flow  at  a  fixed  axial  location  downstream  of  the  nozzle  exit  plane.  A  sche¬ 
matic  diagram  of  the  jet  nozzle  and  the  location  code  is  given  in  Fig.  5.1. 

The  microprocessor  was  programmed  to  obtain  a  fixed  number  of  velocity  measure¬ 
ments  at  each  location  and  then  to  move  automatically  to  the  next  radial  location. 
In  the  data  presented  here,  the  measurement  locations  were  one  millimeter  apart. 
Velocity  distributions  were  obtained  at  axial  positions  spaced  one  centimeter 
apart.  Typical  data  are  presented  in  Figs.  5.2  through  5.20.  The  solid  line 
Is  a  series  of  straight  line  segments  drawn  through  the  mean  velocity  at  each 
radial  point.  The  bars  represent  the  turbulent  velocity  excursions  about 
the  mean  velocity.  For  example,  at  an  axial  position  1  cm  downstream  from 
the  nozzle  exit  the  turbulent  intensity  in  the  potential  core  is  less  than 
2  percent,  but  in  the  mixing  layers  on  the  edges  of  the  exhaust,  the  turbulent 
intensities  are  between  60  and  70  percent. 

As  can  be  seen  in  Fig.  5.2,  the  mean  velocity  Is  relatively  constant 
at  about  65  m/sec.  in  the  core  of  the  jet.  At  axial  locations  further  down¬ 
stream  from  the  exit  plane,  the  mixing  layer  thickness  at  the  edge  of  the  jet 
grows  until  it  finally  reaches  the  jet  centerline.  This  determines  the  length 
of  the  so-called  potential  core  region  which  appeared  to  be  about  13  to  15  cm 
for  the  jet  flow  conditions  presented  here. 

The  distributions  of  turbulent  intensity  at  three  axial  locations 
(nozzle  exit  and  5  and  10  cm  downstream)  are  presented  in  Fig.  5,21.  The 
turbulent  intensities  become  quite  large  in  the  edges  of  the  turbulent  mixing 
layer.  Further  out,  the  intensities  would  of  course  become  very  small. 


The  covariance  (defined  as  (v_,  -  vx)  (vv  -  vy)  /  (N-l)  is  presented  in 
Fig.  5.22  for  three  axial  locations  downstream  from  the  nozzle  exit  plane.  The 
covariance  is  related  to  the  turbulent  shear  stress,  so  that  quite  narrow  zones 
of  intense  shearing  stress  are  present.  For  example  at  the  nozzle  exit,  the 
zones  are  located  at  t  15  mm  from  the  jet  centerline.  The  thickness  of  the  zone 
at  the  nozzle  exit  is  only  3  to  4  mm.  At  an  axial  location  of  5  cm,  there 
appear  to  be  two  zones  of  intense  shearing  stress,  one  edge  at  -  14  mm  and  the 
other  at  +17  mm  and  -18  mm.  There  also  appear  to  be  two  zones  at  an  axial 
location  of  10  cm  downstream  one  at  13  mm  and  15  mm  on  the  positive  side  and 
one  at  -12  mm  and  -16  mm  on  the  negative  side  of  the  centerline. 


Fig.  5.1  Schematic  diagram  of  jet  exhau 
measurement  locations. 


Fig.  5.3  Velocity  distribution  across  subsonic  ]t't. 
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Fig.  5.6  Velocity  distribution  across  subsonic  jut. 
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Fig.  5.7  Velocity  distribution  across  subsonic  jet. 
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Fig.  5.8  Velocity  distribution  across  subsonic  jet. 
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Fig.  5.10  Velocity  distribution  across  subsonic  id. 
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Fig.  5.11  Vi'loeitv  distribution  across  subsonic  let. 
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Fi^.  5.14  Velocity  distribution  across  subsonic  let. 
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Fig.  5.16  Velocity  distribution  across  subsonic  jet 
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5.2  Laser  Veloclmeter  Measurements  in  the  Wake  of  a_  Right  Circular  Cy_l  i  rider 

The  UTSI  Laser  Veloclmeter  System  was  utilized  to  obtain  the  velocity 
distribution  behind  a  0.75  inch  diameter  right  circular  cylinder  In  the 
UTSI  Subsonic  Free  Jet  Wind  Tunnel.  A  schematic  diagram  of  the  free  jet 
nozzle  and  the  cylinder  (viewed  end  on)  is  given  in  Fig.  5.23.  Indicated 
on  the  diagram  is  the  location  code  for  the  laser  probe  volume  which  is 
automatically  scanned  at  a  fixed  axial  location  by  the  LDV  microprocessor 
system. 


A  typical  jet  free  stream  velocity  profile  is  given  in  Fig.  5.23 
for  a  free  stream  velocity  near  65  m/sec.  (which  corresponded  to  a  Mach 
number  of  approximately  0.2). 

As  can  be  seen  in  Fig.  5.25,  the  velocity  at  a  point  downstream  of 
the  cylinder  is  highly  turbulent.  The  data  presented  in  Fig.  5.25  were 
obtained  at  Axial  Position  6,  Scanner  Position  1.  A  statistical  model  to 
treat  these  velocity  spectral  data  is  currently  under  development  at  UTST. 

A  typical  microprocessor  output  listing  for  the  turbulent  velocity 
data  shown  In  Fig.  5.25  is  presented  in  Fig.  5.26.  .Vote  from  the  data 
listing  that  the  mean  velocity  is  10.76  m/sec.  with  a  turbulent  intensity 
of  96  percent . 

As  shown  in  Fig.  5.27,  which  is  the  complete  output  listing  at  a 
typical  point,  the  microprocessor  calculates  the  velocity  and  flow  angle 
for  each  particle  event  and  then  determines  the  mean  flow  velocity  and 
mean  flow  direction  over  the  total  measurement  time  at  each  scan  location. 
Typical  data  for  the  "XV"  scan  mode  are  presented  in  Figs.  5.28  to  5.31. 
Unfortunately,  the  data  for  the  Z-component  ("XZ."  mode)  were  poor  due  to 
spatial  resolution  problems  with  the  third  (or  component.  The  curves 

labeled  "Vv"  represent  flow  in  a  direction  parallel  to  the  cylinder  axis. 

The  data  presented  in  Figs.  5.28  to  5.31  indicate  the  presence  of 
a  reversed  flow  wake  region  behind  the  cylinder.  It  should  he  pointed  out 
that  the  decreases  in  velocity  for  positions  greater  than  0.7  are  due  to 
interaction  with  the  free  jet  mixing  layer.  The  velocity  point  data  can  be 
processed  to  vield  a  velocity  contour  map  as  shown  in  Fig.  5.32.  The  location 
of  the  negative  7.5  m/sec.  contour  is  somewhat  unstable  due  to  the  highly 
turbulent  nature  of  the  flow  behind  the  cylinder. 

The  turbulent  intensities  cx  and  rv  were  calculated  hv  the  micro¬ 
processor  at  each  point  in  the  flow.  Typical  data  for 

cr2x  -  (u  -  if)'  /  XT  ’  are  given  in  Fig.  3.  13. 

A  contour  map  of  rtx  |  s  presented  in  Fig.  5.  1-3.  Note  that  the  highest 
turbulent  intensities  are  generated  near  a  3'  value  of  0.6  inches.  Tt 
is  not  clear  whether  this  region  of  high  turbulence  is  dm-  to  Interaction 
between  the  free  jet  turbulent  mixing  (aver  and  the  cylinder  wake. 
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g.3.2  5  Schematic  diagram  of  jet  exhaust 
measurement  locations. 
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Fig.  5.26  Typical  microprocessor  output  listing. 
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Fig.  5.27  Typical  microprocessor  output  of  1 wo-componcnt  vcTocimetcr . 
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A  t vpical  microprocessor  output  i  ••  presented  in  I  i  >•  .  t  >.  Not  . 
th.it  the  VRMS  is  t  lie  mein  f  I  ow  vole,  itv  in  the  plan.  o!  t  h.  nonet  whi.  h 
is  obtained  from  the  vector  sum  of  Vx  and  \\  which  ore  tin  vi  le,  it  . 

ponents  parallel  and  perpendicular  to  tin  tree  ;el  .  enter!  1  n.  .  r.  .p.ctivi  1  >.  . 

The  respective  STtlMAs  ire  the  velocitv  !  I  no t uat  ions  ' 

—  2  V  -  2  2  i 

[(V  -  V)  ]*  and  [(V  -  V)  /V  ],  respectively.  For  j 

.  \ainp  !  e  ,  a  t  o  ;  .  s  .  :  ; .  ■  .  i  >  .  .  PM  ,  - , 

the  velocitv  s i gma  is  8. iU  m/sec.,  and  the  percent  velocity  1  hwlnot ion  is  f  . 

19.9  percent . 

The  free  stream  velocitv  protile  obtained  upstream  of  t a.  nose  o! 
the  cone  is  presented  in  Tip.  a.  lit.  Shown  on  ea.  h  mean  velocitv  data  -oiut 
is  the  magnitude  of  the  turbulent  velocitv  t  1  no  t  ua  t  i  •■'■l .  Note  that  tin 
velocity  fluctuations  are  quite  1  ante  tot'  the  or.l«  r  of  <•')- (  .  ;  er-a-it  1  at 
the  edges  of  the  flow  where  the  free  jet  mixes  with  the  ambletit  cat.  sunt 
air,  hut  decrease  to  less  than  two  percent  near  the  free  >«•!  .  at.  r in.  . 

Typical  measured  velocitv  pro!  ilos  are  presented  [  n  'in-,.  ..  •' 

through  5.41.  The  procedure  used  to  obtain  the  data  was  to  ,  cat.  t  . 
laser  probe  volume  at  a  fixed  value  ot  X  and  V  (horizontal  ind  v<  rti.  . 
distance)  from  the  cone  nose  tip  iviercmv  point  and  then  the  mi.  topi.  ,  t 
automatically  scanned  the  laser  probe  volume  in  the  .'-direction  i .  at 
the  plane  of  the  paper).  As  noted  in  the  microprocessor  Inpit  .lit  s  -.how: 
in  Fig.  5.35,  the  L.DV  svstem  made  500  particle  v«  locitv  me  isar.  it:  -  it 
each  location  and  then  automatically  moved  on  to  the  tuxi  .ti  n. 

Typical  data  acquisition  times  were  7  to  H  secon.ls  it  each  lo.  .it  i  >u . 

The  vector  velocitv  distribution  in  the  plane  of  t h.  .on.  .nt.: 
line  is  presented  in  Fig.  5.42.  The  velocitv  vector  magnitude  i>  r.  nr.  .  : 
hv  the  length  of  the  vector  while  the  direct  ion  in  the  XV  plan,  i  iv."  1 
the  arrow.  The  axial  velocitv  component  variation  in  the  a  ntrrl  in.  plan, 
is  presented  in  Fig.  5.4).  The  presence  ot  a  vertex  strn.  tnr.  ’a  hind.  •'  . 
cone  can  be  observed  in  Fig.  5.42. 

Due  to  a  programming  error  the  m  i  croprocessor  scantn  .1  tin  '  !■  a  r 
probe  volume  at  the  eompl iment arv  angle  rather  than  scanning  cnrtil.  i  t .  •  1 

the  hodv  surface. 
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Supersonic  Jet  Exhaust  Plume  Measurements 


5.4 


The  UTSI  Laser  Velocimeter  System  operating  in  the  forward  scatter 
mode  (shown  in  the  photograph  in  Fig.  5.44)  was  used  to  measure  the  velocity 
distribution  in  the  exhaust  from  a  Mach  2  nozzle  having  an  exit  diameter  of 
33  mm  (1.3  in.).  During  the  measurements  described  here,  the  plenum  chamber 
pressure  was  111  psia.  Due  to  a  problem  with  the  air  supply  system,  an 
appreciable  amount  of  water  vapor  was  present  in  the  flow.  In  addition,  the 
air  temperature  was  quite  low  (267°K). 

The  velocity  distributions  obtained  for  the  supersonic  nozzle  are 
presented  in  Figs.  5.45  to  5.48  for  X/D  locations  of  0.12,  2.0,  5.0  and  10.0. 

The  velocity  distribution  along  the  jet  centerline  is  presented  in 
Fig.  5.49.  As  expected,  there  are  several  axial  locations  where  the  velocity 
decreases  quite  rapidly  across  shock  waves.  The  data  seem  to  indicate  that 
the  particle  velocities  do  not  go  subsonic  behind  the  shock  waves.  It  is 
thought  that  the  probable  presence  of  water  droplets  and  large  diameter 
particles,  along  with  their  correspondingly  long  velocity  relaxation  distance, 
prevent  equilibration  between  the  gas  and  particle  velocities. 
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Fig.  S.  -j"  Velocity  distribution  along  Mach  2  jet  centerline. 
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5.5  Spectral  Analysis  of  Turbulence 

Comprehensive  theories  tiave  been  developed  for  the  spectral  anal  vs  is 
of  regularly  sampled  (equally  spaced  in  time)  signals.  However,  the  laser 
velocimeter  data  occur  at  random  times  due  to  the  random  arrivals  of 
scattering  particles  in  the  laser  probe  volume. 

A  study  of  the  existing  spectrum  analysis  techniques  for  randomly 
sampled  data  lias  indicated  that  most  are  unsuitable  for  real  time  appl i cat  ions 
and  require  excessive  amounts  of  high-speed  computer  memory. 

During  tin.  previous  US  Army  Research  Office  sponsored  oro',ect,  a 
new  approach  to  spectrum  estimation  from  randomly  sampled  data  records  lias 
been  developed.  The  new  algorithm  is  based  upon  the  use  of  the  method 
Fourier  Series,  the  method  of  least  squares,  and  random  sampling  theory. 

A  comparative  study  indicates  that  it  is  much  faster  than  currently  avail  ah:, 
algorithms  and  can  he  implemented  on  a  small  microprocessor  computer  for  usi 
in  the  semi-real  time  analysis  of  turbulent  velocity  data. 

An  experimental  study  has  been  performed  with  the  result  that  a 
reliable  spectrum  may  be  estimated  with  an  under  sampling  factor  >.f  less 
than  4  (i.e.,  mean  sampling  rate  1/4  of  the  Nvquist  rate). 

Typical  results  obtained  from  computed  simulated  signals  and.  laser 
velocimeter  data  are  presented  in  this  section. 

Figure  5.50  (a)  and  5. 50  (h)  show  -  t  he  lot.-'  a=-,  :  1  and:  a-  s 

characteristics  of  the  algorithm.  A  briok-wal  1  filter,  d  !  t  .•  ;r  j«e  i  >• 
generated  using  the  computer  and  is  used  n  lata  :  •  .  a  I v  >ri t  h-  . 

Computer  simulated  white  noise  did  not  hav.  -'la*  •  .  ■■■■■■•, 

Figure  5.51  and  5.52  depict  no  ise  r.  .  t  i  r  '  :  t  .  s  : :  ;  1  •  t  e 

the  turbulent  velocity  data,  white  noise  i-  add.  '  :  •  a  >•■•  "u 

sine  wave  of  10  Hz.  Figure  5.51  represents  t  h>  t  '  —  tin  •••'  'da  net  s< 

signal.  Figure  5.52  is  the  spectrum  obtained.  in.-  t  I.  a  1  •••rit  lim. 

A  calibration  experiment  was  performed  using  I.'."'  .  u:  is  shewn  in 

Fig.  5.53.  In  this  experiment,  a  loud  speaker  is  1.  •  •  near  the  laser  pr.'b. 
volume  and  is  excited  with  an  acoustic  wave  having  a  frecm-n. v  of  apnroxi- 
matelv  100  Hz  from  a  signal  generator.  The  discrete  velocity  readings  in, 
domain  are  plotted  in  Fig.  5.54.  Figure  5.55  is  the  or  respond  i  ng  s:>e<  t  rum 
from  the  algorithm. 

Considerably  more  work  remains  to  be  done,  hut  the  technique  mav 
permit  the  velocity  power  spectrum  to  he  used  to  determine  the  frequencies 
present  in  the  separated  flow  behind  bodies.  It  is  hoped  that  the  prooediiri 
will  yield  the  vortex  shedding  frequencies  and  the  resulting  turbulence  scale 
1 engths . 
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6.0  CONCLUDING  REMARKS  AND  RECOMMENDATIONS 


Due  to  the  complexity  of  the  laser  velocimeter  system  ultimately 
required  to  obtain  the  required  measurements,  considerably  more  time  and 
effort  than  anticipated  were  required  to  make  the  system  fully  operational. 

In  addition,  flow  quality  problems  associated  with  using  the  NASA  Marshall 
7-inch  Wind  Tunnel  and  the  lack  of  funds  to  support  tests  in  the  AKDC 
Wind  Tunnel  facilities  resulted  in  the  nonavailability  of  a  suitable 
supersonic  test  facility.  Consequently  all  the  planned  tests  were  conducted 
in  a  UTSI  Free  Jet  Facility  which  greatly  restricted  the  size  of  the  model 
which  could  be  employed.  As  a  result,  the  planned  application  to  the 
measurement  of  the  flow-field  generated  by  missile  bodies  at  hi gh-angies- 
of-attack  was  only  partially  completed. 

Some  excellent  experimental  data  were  obtained  for  subsonic  (M  -  (1.2) 
flow  over  right  circular  cylinders  and  for  cones  at  one  high-angl e-of-at tack 
(75°).  While  these  data  are  qualitatively  good,  they  provide  only  a  limited 
quantitative  description  of  the  separated  flow-field  behind  missile  e one- 
cylinder  hodv  combinations.  Considerably  more  data  are  required  to 
characterize  the  nose  vortex  formation  and  separation  process  as  well  as 
the  formation  and  separation  of  body  vortices. 

Now  that  the  UTST  Four-Component  Laser  Velocimeter  and  Microprocessor 
Data  Acquisition  and  Reduction  Systems  are  operational,  it  is  a  straight¬ 
forward  and  cost-effective  procedure  to  obtain  experimental  turbulence 
and  velocity  distribution  data  in  the  body  boundary  layer  and  separated 
wake  regions. 

It  is  recommended  that  an  experimental  measurements  prog. ram  he 
initiated,  using  the  UTST  Laser  Velocimeter  Systems,  to  obtain  data  pertinent 
to  tiie  characterization  of  missile  nose  vortex  formation,  flow  separation 
from  the  missile  body,  and  vortex  growth  and/or  dissipation  in  tile  downstream 
wake  region.  These  data  can  contribute  to  a  better  fundamental  understanding 
of  turbulent  phenomena. 

It  is  further  recommended  that  .in  experimental  program  be  initiated 
to  map  in  a  more  complete  manner  the  wake  region  behind  rone,  cylinder, 
and  cone-cylinder  missile-like  bodies  at  both  subsonic  and  supersonic 
Mach  numbers.  These  data  can  he  usi d  t o  validate  and  improve  current 
missile  high-angle-of-attack  missile  aerodynamics  analytical  models. 
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